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Abstract 

Purpose: The aim of this study is to compare the dose in various soft tissues in brachytherapy with photon emitting 
sources. 

Material and methods: 103 Pd, 125 I, 169 Yb, 192 Ir brachytherapy sources were simulated with MCNPX Monte Carlo 
code, and their dose rate constant and radial dose function were compared with the published data. A spherical phan- 
tom with 50 cm radius was simulated and the dose at various radial distances in adipose tissue, breast tissue, 4-com- 
ponent soft tissue, brain (grey/white matter), muscle (skeletal), lung tissue, blood (whole), 9-component soft tissue, 
and water were calculated. The absolute dose and relative dose difference with respect to 9-component soft tissue was 
obtained for various materials, sources, and distances. 

Results: There was good agreement between the dosimetric parameters of the sources and the published data. 
Adipose tissue, breast tissue, 4-component soft tissue, and water showed the greatest difference in dose relative to the 
dose to the 9-component soft tissue. The other soft tissues showed lower dose differences. The dose difference was also 
higher for 103 Pd source than for 125 1, 169 Yb, and 192 Ir sources. Furthermore, greater distances from the source had higher 
relative dose differences and the effect can be justified due to the change in photon spectrum (softening or hardening) 
as photons traverse the phantom material. 

Conclusions: The ignorance of soft tissue characteristics (density, composition, etc.) by treatment planning systems 
incorporates a significant error in dose delivery to the patient in brachytherapy with photon sources. The error depends 
on the type of soft tissue, brachytherapy source, as well as the distance from the source. 
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Purpose 

Cancer is a main cause of mortality, and the number of 
cancer diagnoses annually is more than ten million world- 
wide [1]. Among the currents cancer treatment modalities, 
the most common are chemotherapy, radiotherapy, and 
surgery [2]. As a requirement for a successful radiothe- 
rapy, it is essential to minimize the discrepancies between 
calculated dose distributions and delivered doses. To this 
end, the accurate calculation of dose distributions in the 
treatment planning process is paramount [3,4]. Tissue and 
water equivalent materials are routinely used in quality 
assurance and calibration processes involved in diagnostic 
radiology and radiotherapy. However, from a radiation 
physics point of view, currently available tissue equiva- 
lent materials show limitations in mimicking characteris- 
tics of the real tissues in low and high energy ranges [5]. 



Treatment planning systems (TPSs), which are used in 
radiotherapy, utilize correlations based on in-water mea- 
surements for evaluation of doses in the volume of interest 
and in other organs. TPSs exhibit errors in estimation of or- 
gan doses in some cases. For example, TPSs overestimate 
the dose inside heterogeneity and at points beyond it. It 
has been shown that the algorithms, which are employed 
for computations in the TPSs are not able to calculate 
the dose variation in heterogeneous regions as accurately 
as Monte Carlo methods [6-9]. In a number of previous 
studies, the dose distribution in various tissue-equivalent 
materials was calculated. Yazdani and Mowlavi have de- 
termined dosimetric parameters of a 131 Cs brachytherapy 
source (model CS-1 Rev.2), based on a Monte Carlo study. 
They have also calculated relative dose distributions in 
water and soft tissue phantoms. It was concluded that the 
dose deposition from photons in high gradient regions 
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can be calculated accurately by Monte Carlo techniques. 
The results of the study can be used as input for TPSs 
using this source [10]. Mowlavi and Yazdani performed 
a similar study on a high dose rate 169 Yb source (model 
X1267), yielding similar conclusions [11]. Luxton, based 
on the EGS4 Monte Carlo code, has developed a mod- 
el for calculation of dose rate in water for a low-energy 
brachytherapy source from the measurement data of dose 
rate to water within a solid water phantom. The model 
was used for calculation of the dose rate at various dis- 
tances from a source in a thin shell of water within the sol- 
id phantom. It was concluded that by comparing the dose 
distribution with that calculated for a homogeneous water 
phantom, correction factors for determination of dose rate 
in homogeneous water medium from the measurements 
within the solid phantom can be derived. The correction 
factors were calculated for acrylic, solid water (WT1), and 
RW-1 phantoms with photon spectra of a 103 Pd and two 
125 I sources used as input. Among the three phantom ma- 
terials evaluated for these three sources, measurements in 
RW-1 phantom were most similar to the corresponding 
values in a water phantom [12]. Furthermore, there are 
also other studies on various tissue-equivalent materials 
when phantoms irradiated by brachytherapy sources to 
evaluate their water equivalence and their dose distribu- 
tion differences [13,14]. 

Since the presence of inhomogeneities in a phantom 
affect the dose distribution, some previous studies have 
been focused on dose determination in the presence of in- 
homogeneity. Chandola et al. performed a study to find 
out the dose differences in the presence of inhomoge- 
neities, in a water phantom, using a high dose rate 192 Ir 
brachytherapy source. The source, a water phantom, 
and inhomogeneities were simulated using Monte Car- 
lo simulation EGSnrc code. The relative dose difference 
along the transverse axis of the source was 5.5-6.5% high- 
er and 4.5-5% lower in the presence of air and cortical 
bone, respectively. The results were in good agreement 
with other data in the literature [15]. Cazeca et al. used 
MCNP5 code to calculate the dose rate within planning 
target volume of a MammoSite balloon® (RTS Cytyc 
Corp, Marlborough, MA, USA) dose delivery system for 
evaluation of the effects of breast-air and breast-lung in- 
terfaces on the dose within this volume. An average fe- 
male chest phantom and a semi-infinite water phantom, 
and 169 Yb and 192 Ir high dose rate brachytherapy sources 
were simulated. The ratio of the calculated dose rate in 
the planning target volume for the chest and the water 
phantom in the point closest to the breast-air interface 
was less than that in the point closest to the breast-lung 
interface by 11.4% and 4% for the 169 Yb, and 192 Ir sources, 
respectively. The results imply that to avoid error in dose 
delivery into the patient, the software which is used for 
dose rate calculation in MammoSite balloon® treatments 
should account for patient anatomy and the density of 
the surrounding materials in its dosimetric analysis [16]. 
In a Monte Carlo study by Gialousis et al., the dose dif- 
ference originating from tissue inhomogeneity was calcu- 
lated. The results showed that a relative underestimation 
in dose by TPS calculations existed in tissues close to the 



radioactive source. The existence of lung tissue instead of 
normal tissue resulted in relative dose increase of 8% at 
4-cm from the source. Furthermore, the relative increase 
was 2.1% for the lung and amounted to 6.8% and 7.6% 
for the organs at risk such as heart and bone marrow, re- 
spectively [17]. Hsu et al. has reported dose distributions 
around a 192 Ir brachytherapy source in various phantoms 
using three Monte Carlo codes. Furthermore, measure- 
ments were performed for verification of the Monte Carlo 
results. The results have shown that radial dose functions 
were influenced in bone tissue. Dose differences ranging 
from 0.6% to 14.4% were observed between the values 
in the homogeneous solid water phantoms and the solid 
water, which included lung interfaces. The dose differ- 
ence between the bone phantoms and bone-lung interfac- 
es ranged from 4.1% to 15.7%. These results imply that 
there are dose distribution differences in water, bone, 
lung, and inhomogeneous phantoms. It was concluded 
that clinical parameters do not provide sufficient dose 
calculation accuracy for different materials, and to im- 
prove the brachytherapy treatment quality calculations, 
TPSs should incorporate material density [18]. 

In a study by White et al., the dosimetric impact of 
trace elements in human tissues was determined for low 
energy photon brachytherapy sources. White et al. have 
used Monte Carlo calculations to investigate the effect 
of trace elements in normal and cancerous tissues on the 
dose in brachytherapy with low-energy photon sources. 
Their results indicate that the presence of trace elements 
in tissue result in dose differences, which depend on the 
atomic number and concentration of the trace elements. 
While trace elements of low atomic number have a negli- 
gible effect in all tissues, higher atomic number elements 
had larger effects, which are higher than 3%. Further- 
more, the dose distributions in cancerous and healthy 
prostate tissues and trace element-free composition 
were significantly different. It was mentioned that with 
low-energy sources, due to the non-negligible effect on 
the dose in tissues, the trace elements are taken into ac- 
count as a source of uncertainty in dose calculations and 
further investigations are needed for accurate determina- 
tion of the trace composition effect in various soft tissues 
[19]. While various studies have evaluated the dose in 
various tissue-equivalent material phantoms and in one 
case the effect of in-tissue trace elements was evaluated, 
to the best of our knowledge the difference between the 
relative and absolute dose in various soft tissues has not 
been evaluated, from a comparison point of view. 

The aim of this study is to compare the dose in wa- 
ter (as a dosimetric material reference) and various soft 
tissues with ICRU 9-component soft tissue [20] in bra- 
chytherapy with photon emitting sources. 

Material and methods 

Brachytherapy sources 

Four brachytherapy sources were simulated in this 
study: 103 Pd (Prospera Med3633 ideal, Brachytherapy 
Services, Inc. CA, USA), 125 I (Amersham QncoSeed 6702, 
Little Chalfont, United Kingdom), 169 Yb (4140, Implant 
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Sciences Corporation, Wakefield, MA, USA), 192 Ir (Gam- 
maMed 12i, Varian Medical Systems, Inc. Palo Alto, CA, 
USA) brachytherapy sources were simulated with MCN- 
PX Monte Carlo code (version 2.4.0) [21]. The Med3633 
103 Pd source is composed of two polystyrene spheres, each 
having 0.560 mm diameter coated with a negligible thick- 
ness of radioactive 103 Pd. The polystyrene spheres are po- 
sitioned on either side of two gold/ copper alloy spheres 
(markers) with 0.560 mm diameter. The seed has a titani- 
um encapsulation, which is in the form of a cylinder with 
outer and inner diameter of 0.810 mm and 0.710 mm, re- 
spectively. The average thickness of the weld is 0.100 mm. 
The centers of the source spheres are located at ± 1.807 mm 
and ± 1.084 mm. The centers of the markers are located 
at ± 0.361 mm distances from the center of the seed. The 
source has overall length of 4.70 mm, while its active 
length is 4.20 mm [22]. 

The Amersham OncoSeed 6702 125 I source includes 
three resin spheres with 1.2 g/cm 3 density and molecular 
formula of C 12 H lg NCl. The diameter of the resin spheres 
is 0.600 mm. A negligible thickness of radioactive 125 I is 
coated on the spheres. The encapsulation of the source 
is composed of a titanium tube with 0.050 mm thickness 
on its walls having an outer diameter of 0.800 mm. There 
are two end welds with 0.500 mm thickness, which are in 
the form of 0.400 mm hemispheres located on top of solid 
cylinders with 0.400 mm radius and 0.100 mm thickness. 
The active length of the source is 3.30 mm and its overall 
length is 4.50 mm [23]. 

The HDR 4140 169 Yb source consists of a ytterbium 
oxide core (6.9 g/cm 3 ) with length of 3.60 mm and dia- 
meter of 0.73 mm. The active core has a stainless steel 
capsule (AISI 306 with 7.80 g/cm 3 density). The end weld 
consists of a hemispherical part with 0.45 mm radius and 
a solid cylindrical part with 0.050 mm thickness. A hol- 
low cylindrical section with 4.30 mm length is attached to 
the end weld. The inner and outer diameters of the hol- 
low cylindrical section are 0.73 mm and 0.90 mm, respec- 



tively. This hollow portion is connected to a cylindrical 
section with 0.60 mm thickness and 0.90 mm diameter. 
The source's cable is in the form of a solid cylinder made 
of stainless steel (AISI 306 with density of 6.90 g/cm 3 ) 
with diameter of 0.90 mm and length of 1.95 mm [24]. 

The GammaMed 12i 192 Ir source consists of an ac- 
tive 192 Ir core and a stainless steel capsule. The length of 
the active core is 3.50 mm and its diameter is 0.70 mm. 
The capsule has 1.10 mm diameter and is made of AISI 
31 6L stainless steel with a density of 7.8 g/cm 3 . The tip 
of the source, which is part of its encapsulation, has the 
form of a cone with 0.143 mm height and opening angle 
of 75°. This section is attached to solid cylindrical section 
with a length of 0.717 mm followed by a hollow section 
with length and inner diameter of 3.60 mm and 0.70 mm, 
respectively. Attached to the hollow section is a 0.5 mm 
long solid cylinder. A 6 cm stainless steel cable (AISI 304 
with density of 5.6 g/cm 3 ) is connected to the source [25]. 

A schematic diagram of the geometry of the Pros- 
pera Med3633 103 Pd, Amersham OncoSeed 6702 125 1, 4140 
169 Yb and GammaMed 12i 192 Ir brachytherapy sources 
is illustrated in Figure 1. Photon energy spectrum of 
103 Pd, 125 I, 169 Yb and 192 Ir is illustrated in Table 1 [26-28]. 
Some of these radionuclides, such as 169 Yb, are not as wide- 
ly used as the 192 Ir source. However, the four brachythe- 
rapy sources evaluated in this study, spanning a relatively 
wide energy range, were selected to study the energy de- 
pendency of the results. 

Dosimetric parameters of sources 

Dose rate constant and radial dose function of the 
103 Pd, 125 I, 169 Yb, and 192 Ir brachytherapy sources were 
calculated based on updated AAPM task group No. 43 
(TG-43U1) formalism [26]. For calculation of air kerma 
strength, a number of torus cells containing dry air were 
defined at distances ranging from 1-50 cm from each 
source and air kerma rate in the cells was calculated using 
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Fig. 1. A schematic diagram illustrating the geometry of (A) Prospera Med3633 103 Pd, (B) Amersham OncoSeed 6702 125 I, 
(C) 4140 169 Yb, and (D) GammaMed 12i 192 Ir brachytherapy sources used in this study. The dimensions are not to a real scale 
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Table 1. Photon energy spectrum for 125 


1 [26], 103 Pd [26], 169 Yb [27], and 


192 lr [28] radionuclides 




125| 


10; 


; Pd 


169 Y b 


192 


Ir 


Energy (keV) Intensity (%) 


Energy (keV) 


Intensity (%) 


Energy (keV) 


Intensity (%) 


Energy (keV) 


Intensity (%) 


27.202 0.405 


22.074 


0.224 


49.77 


53.2 


61.49 


1.20 


27.472 0.756 


20.216 


0.423 


50.74 


94.0 


63.00 


2.07 


30.980 0.201 


22.717 


0.104 


57.60 


29.5 


65.12 


2.65 


31.877 0.0438 


23.312 


0.0194 


59.10 


8.2 


66.83 


4.53 


35.4919 0.0466 


39.755 


0.000683 


63.10 


44.2 


71.08 


0.24 




62.51 


0.0000104 


93.62 


2.6 


71.41 


0.46 




294.52 


0.000028 


109.78 


17.5 


73.36 


0.16 




357.46 


0.000221 


118.19 


1.9 


75.37 


0.53 




497.054 


0.0000401 


130.52 


11.3 


75.75 


1.03 








177.21 


22.2 


77.83 


0.37 








197.96 


35.8 


136.34 


0.18 








261.08 


1.7 


201.31 


0.47 








307.74 


10.1 


205.80 


3.30 












283.27 


0.26 












295.96 


28.67 












308.46 


30.00 












316.51 


82.81 












374.49 


0.72 












416.47 


0.66 












468.07 


47.83 












484.58 


3.18 












489.04 


0.44 












588.58 


4.52 












604.41 


8.23 












612.47 


5.31 












884.54 


0.29 



Mean energy = 28.37 keV Mean energy = 20.74 keV Mean energy = 92.7 keV Mean energy = 350 keV 



F6 tally. Air kerma rate x r 2 was plotted versus r (cm) 
and its average in the relatively flat region of the plot (2- 
30 cm) was used toward calculation of air kerma strength. 
Air torus cells were located in a vacuum medium. 
The Monte Carlo calculation type A statistical uncertain- 
ties were less than 0.45%. Type A statistical uncertainty 
refers to random error [29] . Dose rate in water was scored 
using *F8 tally, and then dose at 1 cm distance from the 
source was divided to the value of air kerma strength and 
dose rate constant (cGy h" 1 U" 1 ) was calculated. The water 
phantom in this calculation was a sphere with 50 cm radi- 
us. The Monte Carlo errors in the tally cell for calculation 
of the dose were maximally 0.32%. 

Geometry function was calculated using the line- 
source approximation and then was used for calculation 
of radial dose function. For the calculation of radial dose 
function, a water phantom as a sphere with 50 cm radi- 
us was defined and dose values at distances ranging from 
0.5-7 cm was calculated by utilizing *F8 tally. The calcula- 
tion of radial dose function was based on equation No. 6 



of TG-43U1 formalism [26]. In this section, for all the 
sources, 1.5 x 10 8 photons were scored. The maximum 
Monte Carlo error in the tally cells was 4.22% for the case 
of calculation of radial dose function. 

The effect of tissue composition 

Mass density and composition (fraction by weight) 
of adipose tissue, breast tissue, 4-component soft tissue, 
brain (grey/white matter), muscle (skeletal), lung tissue, 
blood (whole), and 9-component soft tissue were adopted 
from the report No. 44 of the International Commission 
on Radiation Units and Measurements (ICRU) [20] . Mass 
density and the composition of the mentioned tissues and 
water are listed in Table 2. Each source and phantom ma- 
terial was simulated using MCNPX (version 2.4.0) code 
[21]. To have an acceptable level of Monte Carlo type A 
statistical uncertainty within a reasonable computation 
time, the absorbed dose was determined by calculation 
of the photon flux and then multiplication of the flux 
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Table 2. Mass density and composition (fraction by weight) of various tissues evaluated in this study. The data 
are extracted from ICRU-44 report 





Adipose 
tissue 


Breast 
tissue 


Soft tissue 
(4-compo- 
nent) 


Brain 
(grey/ 
white 
matter) 


Muscle 
(sKeietaij 


Lung 
tissue 


Blood 
(whole) 


Soft tissue 
(9-compo- 
nent) 


Water 


Density (g/cm 3 ) 0.95 


1.02 


1.00 


1.04 


1.05 


1.05 


1.06 


1.06 


1.00 


H 


0.114000 


0.106000 


0.101174 


0.107000 


0.102000 


0.103000 


0.102000 


0.102000 


0.111111 


C 


0.598000 


0.332000 


0.111000 


0.145000 


0.143000 


0.105000 


0.110000 


0.143000 




N 


0.007000 


0.030000 


0.026000 


0.022000 


0.034000 


0.031000 


0.033000 


0.034000 




0 


0.278000 


0.527000 


0.761826 


0.712000 


0.710000 


0.749000 


0.745000 


0.708000 


0.888889 


Na 


0.001000 


0.001000 




0.002000 


0.001000 


0.002000 


0.001000 


0.002000 




P 




0.001000 




0.004000 


0.002000 


0.002000 


0.001000 


0.003000 




S 


0.001000 


0.002000 




0.002000 


0.003000 


0.003000 


0.002000 


0.003000 




CI 


0.001000 


0.001000 




0.003000 


0.001000 


0.003000 


0.003000 


0.002000 




K 








0.003000 


0.004000 


0.002000 


0.002000 


0.003000 




Fe 














0.001000 







with mass energy absorption coefficient of the soft tissue. 
For this purpose, the energy flux (MeV/cm 2 ) per source 
particle at various distances from the source in a 50 cm 
radius spherical phantom containing the soft tissue was 
calculated using the *F4 tally, and then was converted to 
the absorbed dose by multiplication with the mass energy 
absorption coefficients for the soft tissue extracted from 
the NIST database [30]. Since the mass energy absorption 
coefficient depends on the photon energy, the energy was 
divided in various bins and the flux multiplication was 
performed at each energy bin. This calculation was re- 
peated for each soft tissue using mass energy absorption 
coefficient of that tissue. Each input file for 103 Pd source 
was run for 3 x 10 8 photon histories and the maximum 
Monte Carlo type A statistical uncertainty was 4.26%. 
The input files for other sources were run for 1.5 x 10 8 



photons, and the maximum Monte Carlo type A statisti- 
cal uncertainty was 1.41%. The number of photons scored 
for 103 Pd source was more than for other sources, because 
for 103 Pd source, the MC errors in the tally cells were be- 
ing reduced with a slower speed. The relative dose in adi- 
pose tissue, breast tissue, 4-component soft tissue, brain 
(grey/white matter), muscle (skeletal), lung tissue, blood 
(whole), and water, to dose in 9-component soft tissue 
was calculated. 

Mass energy absorption coefficient and effective ato- 
mic number (Z eff ) for total photon-energy absorption of 
adipose tissue, breast tissue, 4-component soft tissue, brain 
(grey/white matter), muscle (skeletal), lung tissue, blood 
(whole), 9-component soft tissue, and water at 20, 30, 100, 
and 300 keV photon energies are listed in Table 3. Among 
a number of values in a relatively wide range of photon 



Table 3. Mass energy absorption coefficient values (g/cm 2 ) [30] and Z eff values for total photon-energy absorp- 
tion of various tissues and water [31] in 20, 30, 100, and 300 keV photon energies 



Photon energy 




Adipose 


Breast 


Soft 


Brain 


Muscle 


Lung 


Blood 


Soft 


Water 


(keV) 




tissue 


tissue 


tissue 


(grey/ 


(skele- 


tissue 


(whole) 


tissue 












(4-com- 


white 


tal) 






(9-com- 












ponent) 


matter) 








ponent) 




20 


Pen 


0.32510 


0.43940 


0.50700 


0.57060 


0.56380 


0.57400 


0.58310 


0.56630 


0.55030 


Zeff 


4.9032 


5.4105 


5.7460 


5.8629 


5.9085 


5.9409 


5.9873 


5.9169 


5.7826 


30 


Pen 


0.09495 


0.12600 


0.14380 


0.16290 


0.16100 


0.16350 


0.16690 


0.16160 


0.15570 


Zeff 


4.1280 


4.6653 


5.0403 


5.0971 


5.1532 


5.1749 


5.2157 


5.1590 


5.0282 


100 


Pen 


0.02433 


0.02478 


0.02501 


0.02558 


0.02544 


0.02550 


0.02559 


0.02545 


0.02546 


^eff 


3.1286 


3.3623 


3.5490 


3.4549 


3.5353 


3.5379 


3.5542 


3.5364 


3.4313 


300 


Pen 


0.03194 


0.03173 


0.03161 


0.03178 


0.03164 


0.03167 


0.03164 


0.03164 


0.03192 


Zeff 


3.0779 


3.2918 


3.4648 


3.3627 


3.4422 


3.4432 


3.4569 


3.4429 


3.3432 



Journal of Contemporary Brachytherapy (2014/volume 6/number 1) 



Effect of tissue composition on dose distribution in brachytherapy 



59 



energies, only those at 20, 30, 100, and 300 keV energies 
were presented in Table 3. Selection of these energies 
was based on the average photon energies emitted by 
the 103 Pd (average 20.74 keV (TG43-U1 [26]), 125 I (average 
28.37 keV (TG43U1 [26]), 169 Yb (average 92.7 keV), and 
192 Ir (average 350 keV) radionuclides. 

Calculation of relative and absolute dose 
for the sources 

Relative dose was defined as the ratio of dose in the 
adipose tissue, breast tissue, 4-component soft tissue, 
brain (grey/ white matter), muscle (skeletal), lung tissue, 
blood (whole), and water, to the dose of 9-component soft 
tissue (as a reference tissue in most studies). These values 
were calculated in torus cells at 0.2, 0.3, 0.5, 1, 2, 3, 4, 5, 
7, 10, 12, and 15 cm radial distances relative to the afore- 
mentioned sources to examine how the trend of relative 
dose changes with the distance. 

Relative dose with respect to the dose of 9-component 
soft tissue alone may not have much significance, because 
its value changes with the energy of the source and also 
with radial distance from each source. Furthermore, be- 
cause of a relatively high dose gradient near the source, 
this value alone would not be clinically meaningful. To 
have an interpretation of the variation of the dose rate 
with these parameters, the value of dose rate (cGy h _1 U" 1 ) 
in the adipose tissue, breast tissue, 4-component soft tis- 
sue, brain (grey/ white matter), muscle (skeletal), lung 
tissue, blood (whole), 9-component soft tissue, and water 
was calculated at various radial distances from the sourc- 
es as well. The method for calculation of dose in these tis- 
sues was the same as the aforementioned method. In other 
words, the flux was calculated and then converted to dose 
by multiplication by mass energy absorption coefficient of 
each soft tissue. The Monte Carlo output (MeV/g per pho- 
ton) was converted to dose rate (cGy h _1 U _1 ) as described 
below: 

Dose rate (cGy h" 1 U" 1 ) = MC output (MeV/g per 
photon) x 106 (eV/MeV) x 1.602 x 10' 19 (J/eV) x 1000 
(g/kg) x 1 (dis/s per Bq) x source activity (Bq/per U) x 
x photon yield of the source (photons/ dis) x 3600 (s/h) x 
x 100 (cGy/Gy). (1) 

Calculation of photon spectrum 

In order to interpret the results of dose differences 
in various soft tissues at various radial distances from 
the source, it is useful to know the change in photon spec- 
trum at various distances from the sources. Therefore, 



the photon spectra at 0.2, 1, 5, 10, and 15 cm from the 
103 Pd, 125 I, 169 Yb, and 192 Ir sources were calculated. For 
this purpose, torus cells were defined at the stated dis- 
tances from each source and the photon flux in the cells 
were calculated using F4 tally. The flux was calculated 
in various energy bins with 20 keV energy intervals for 
103 Pd, 169 Yb and 192 Ir sources. In the case of 125 I, energy 
bins were assigned 5 keV. To have an acceptable level of 
Monte Carlo error in the tally cells in calculation of pho- 
ton flux, the importance of photons and electrons in the 
tally cells was set as 50, while the importance in other 
cells was defined equal to 1. The input files were run for 
2 x 10 9 particles for 103 Pd, 125 1, 169 Yb, and 5 x 10 8 for 192 Ir 
source. However, a Monte Carlo error of about 30% was 
observed in flux calculation for the farthest cell for 103 Pd 
source, the error of tally calculation for other cells and 
sources was significantly less than this level. The energy 
spectrum was plotted for the 0.2, 1, 5, 10, and 15 cm dis- 
tances for each brachytherapy source. 

Results 

TG-43 dosimetric parameters 

Dose rate constant and radial dose function values 
for the 103 Pd, 125 1, 169 Yb, and 192 Ir sources obtained in this 
study are listed in Table 4 and Table 5. These tables also 
include the corresponding values from other studies for 
the same source models, and also the percentage differ- 
ences (%) between the data from this study and the lite- 
rature. 

The effect of tissue composition 

The relative dose with respect to 9-component soft 
tissue dose for adipose tissue, breast tissue, 4-component 
soft tissue, brain (grey /white matter), muscle (skeletal), 
lung tissue, blood (whole), and water in the case of 103 Pd, 
125 1, 169 Yb, and 192 Ir brachytherapy source models is pre- 
sented in Table 6. The relative dose here means the ra- 
tio of dose in soft tissue (for example adipose) to dose 
in 9-component soft tissue. The relative dose values were 
listed in this table for various radial distances from each 
source ranging from 0.2-15 cm. 

The dose rate (cGy h" 1 U" 1 ) in various soft tissues for 
the four brachytherapy sources at various radial distanc- 
es from each source is listed in Table 7. While the data in 
the 0.2-4.0 cm distances were rounded to two decimals 
in this Table, to show the differences in a number of data 
points with low dose values, the data in these points were 



Table 4. Dose rate constant (cGy h 1 U : ) for the 103 Pd, 125 l, 169 Yb, 192 lr brachytherapy sources and the percentage 
difference (Diff. [%]) between the values presented here with previously published data 

Source (model) Present study Other study (Reference) Diff. (%) 

103 Pd (Prospera Med3633 ideal) 0.652 0.672 (Rivard et ai [32]) 2.93 

125 l (Amersham OncoSeed 6702) 1.069 1.036 (TG-43U1 [26]) -3.22 

169 Yb(4140) 1.191 1.186 (Taylor et al. [33]) -0.43 

192 lr (GammaMed 12i) 1.139 1.122 (Perez-Calatayud etai [34]) -1.48 
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presented up to four decimals. Figure 2 demonstrates 
the dose rate (cGy h -1 U" 1 ) versus radial distance (cm) for 
various soft tissues and 103 Pd, 125 1, 169 Yb, and 192 Ir source 
models. In this figure, dose rate values were only plotted 
to 2 cm from each source, because dose rate values be- 
yond this distance may not include any important infor- 
mation on the graphs. 

Photon spectrum at various distances front 
the sources 

The photon spectra at 0.2, 1, 5, 10, and 15 cm distanc- 
es from the 103 Pd, 125 I, 169 Yb, and 192 Ir source models are 
illustrated in Figure 3. The spectra plots are the normal- 
ized number of photons per history to the total number 
of photons inside cells defined at mentioned distances 
versus photon energy (keV) in 9-component soft tissue. 

The average photon energies on the spectra for 103 Pd, 
125 I; i69 Y b, and 192 Ir sources at 0.2, 1, 5, 10, and 15 cm 
distances from the sources in a 9-component soft tissue 
phantom are presented in Table 8. The data in this table 
can be used toward interpretation of the variation of rel- 
ative dose in various soft tissues with distance from the 



Discussion 

In the present study, four photon emitting brachythe- 
rapy sources were simulated and validation simulations 
were used to evaluate the effect of tissue composition on 
dose for a number of tissue types. Agreement within 2-5% 
for anisotropy functions were obtained, compared to the 
supplement to the of TG-43 updated report (TG-43U1S1) 
[35]. Considering this criterion as excellent agreement, 
it can be seen in Table 4 that dose rate constant values 
of the present study are in excellent agreement with the 
previously published data. Ignoring some exceptions, 
the calculated radial dose functions for 103 Pd, 125 I, 169 Yb, 
and 192 Ir sources (Table 5) at almost all the data points 
show excellent agreement when compared to the radial 
dose functions reported in other studies. The calculated 
dosimetric parameters were used as the reference values 
for the simulated source models. 

The results of the present study (Table 6 and Table 7) 
indicate that the differences in the composition of various 
soft tissues can influence the dose distribution. The effect 
is considerable and under similar conditions, the dose dis- 
tribution in each soft tissue differs from the others. This 
effect implies that by equating various tissues in dose cal- 
culations with TPSs in brachytherapy by photon emitting 
sources, especially by lower energy sources, a significant 
error is introduced. A similar effect is observed for wa- 
ter as a tissue-equivalent material. The values of relative 
dose in Table 6 and the absolute dose values in Table 7 
for water indicate that the dose in water significantly dif- 
fers from the dose in 9-component soft tissue. This fact 
can be justified by comparing the mass energy absorption 
coefficient, and Z eff of water and 9-component soft tissue 
presented in Table 3. This effect implies that dosimetry 
using water as a tissue equivalent medium, and extension 
of the in-water dose values for soft tissue, incorporates 
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Table 7. Dose rate (cGy h 1 U in various soft tissues for 103 Pd, 125 l, 169 Yb, and 192 lr brachytherapy sources 
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0.0587 




7.0 


0.0279 


0.0321 


0.0292 


0.0297 


0.0296 


0.0297 


0.0296 


0.0296 


0.0297 




10.0 


0.0133 


0.0149 


0.0133 


0.0133 


0.0132 


0.0132 


0.0131 


0.0132 


0.0134 




12.0 


0.0087 


0.0096 


0.0086 


0.0085 


0.0083 


0.0083 


0.0082 


0.0082 


0.0086 




15.0 


0.0050 


0.0054 


0.0048 


0.0046 


0.0045 


0.0045 


0.0044 


0.0045 


0.0047 
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Table 7. Cont. 



r(cm) 


Adipose 
tissue 


Breast 
tissue 


Soft 
tissue 
(4-compo- 
nent) 


Brain 
(grey/ 
white 
matter) 


Muscle 
(skeletal) 


Lung 
tissue 


Blood 
(whole) 


Soft 
tissue 
(9-compo- 
nent) 


Water 










Dose rate (cGy h-1 U-l) 








192 lr 0.2 


28.38 


28.20 


28.13 


28.30 


28.16 


28.20 


28.16 


28.16 


28.42 


0.3 


12.58 


12.50 


12.47 


12.54 


12.48 


12.50 


12.48 


12.48 


12.60 


0.5 


4.53 


4.50 


4.49 


4.52 


4.50 


4.50 


4.50 


4.50 


4.54 


1.0 


1.14 


1.13 


1.13 


1.13 


1.13 


1.13 


1.13 


1.13 


1.14 


2.0 


0.29 


0.29 


0.28 


0.29 


0.28 


0.29 


0.28 


0.28 


0.29 


3.0 


0.13 


0.13 


0.13 


0.13 


0.13 


0.13 


0.13 


0.13 


0.13 


4.0 


0.07 


0.07 


0.07 


0.07 


0.07 


0.07 


0.07 


0.07 


0.07 


5.0 


0.0459 


0.0461 


0.0457 


0.0460 


0.0458 


0.0458 


0.0458 


0.0458 


0.0461 


7.0 


0.0231 


0.0233 


0.0230 


0.0231 


0.0230 


0.0230 


0.0230 


0.0230 


0.0232 


10.0 


0.0136 


0.0138 


0.0135 


0.0136 


0.0135 


0.0135 


0.0135 


0.0135 


0.0137 


12.0 


0.0108 


0.0110 


0.0108 


0.0108 


0.0107 


0.0107 


0.0107 


0.0107 


0.0109 


15.0 


0.0072 


0.0074 


0.0072 


0.0072 


0.0071 


0.0071 


0.0071 


0.0071 


0.0073 



a non-negligible error in dose delivery especially in points 
far from the source model. By introducing organ specific 
density and composition into the TPS and dose calcula- 
tion, it would be possible to avoid this error. Since the re- 
port No. 24 of ICRU has recommended that the accuracy 
required for tumor dose delivery in radiotherapy is ± 5% 
[36] the above effects are considerable, and a TPSs, which 
incorporating the differences in various soft tissues are 
capable of calculating dose distribution more accurately. 

As it can be seen from the data in Table 6 and Table 7, 
while 103 Pd, 125 I, and 169 Yb source models in adipose tis- 
sue, breast tissue, 4-component soft tissue, and water il- 
lustrate greatest difference in dose relative to the dose to 
the 9-component soft tissue, respectively, the other soft 
tissues show lower dose differences. Although the dose 
differences for other soft tissues are lower, this fact is not 
negligible. A detailed discussion on the effect of materi- 
al composition on dose distribution (by considering the 
weight fraction of each element in the material) may be 
difficult, but the composition effect can be incorporated 
by considering the effective atomic number (Z eff ) of the 
soft tissues and water. A comparison of Z e£f and mass 
energy absorption coefficient of water and soft tissues as 
presented in Table 3. This indicates that, with some ex- 
ceptions, the Z eff and mass energy absorption coefficient 
for adipose, breast and 4-component soft tissue are lower 
than other soft tissues. This can explain the higher dose 
differences of these tissues compared to the other tissues. 
However, it should be noted that the data in Table 3 and 
Table 4 are related to only a number of selected pho- 
ton energies, but all the photon energies emitted by the 
source models effect the dose differences for various tis- 
sues (Tables 6 and 7). As it can be seen in Table 2, Z eff and 
mass energy absorption coefficient of soft tissues in high- 
er energies have similar values, for example in 300 keV 



Table 8. Average photon energy for 103 Pd, 125 1, 159 Yb, 
and 192 lr sources at 0.2, 1, 5, 10, and 15 cm distances 
in a 9-component soft tissue phantom 



Source 


Radial distance from the source (cm) 




0.2 


1 


5 


10 


15 


103p d 


38.74 


37.30 


38.55 


75.04 


128.70 


125| 


31.16 


30.79 


30.48 


30.77 


31.26 


169 Yb 


114.31 


108.90 


95.80 


90.77 


88.70 


192, r 


356.50 


330.65 


233.09 


180.23 


154.62 



the maximum difference in Z e£f and mass energy absorp- 
tion coefficient of soft tissues are only 0.3869 and 0.00033, 
respectively. Therefore, values around 1 for relative dose 
with respect to 9-component soft tissue can be predictable 
in the case of 192 Ir source model. 

The results showed the largest relative dose differ- 
ence occurred for low energy source models ( 103 Pd and 
125 I), Table 6. This effect can be justified with considering 
this fact that photoelectric effect has a role in absorption 
of radiation in soft tissue. The probability of this effect 
decreases with increase of photon energy. And since the 
photoelectric effect is inversely proportional to the third 
power of photon energy (1/E 3 ), the absorbed dose dif- 
ference with 9-component soft dose is higher for lower 
photon source models, such as 103 Pd and 125 I, and de- 
creases with energy to its lowest values for 192 Ir source. In 
a study by Melhus and Rivard [37], the effect of phantom 
size, phantom material, mode of radiation transport, vol- 
ume averaging, and Monte Carlo tally type were studied. 
125 I and 103 Pd sources showed the maximum sensitivity to 
phantom material between tissue substitutes and water, 
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Adipose tissue M Breast tissue —A — Soft tissue (4-component) — H — Brain (grey/ white matter) 
«- Muscle (skeletal) — t— Lung tissue — ♦ — Blood (whole) — # — Soft tissue (9-component) Water 

Fig. 2. Dose rate (cGy h" 1 U" 1 ) in various soft tissues for: (A) 103 Pd, (B) 125 I, (C) 169 Yb, and (D) 192 Ir brachytherapy sources 



by a factor of 1.4 and 2.0 at 9 cm distance from the source, 
respectively. Furthermore, it has been shown that there 
are ±5% differences in dose distributions for 137 Cs, 192 Ir, 
and 169 Yb high-energy photon emitting sources between 
water and tissue-equivalent materials at 20 cm distance 
from the source. Due to the differences in methodologies 
between the present study and the study by Melhus and 
Rivard, it may not be appropriate to have a point by point 
comparison between the results of these two studies. 
A difference between these two studies is that in the pres- 
ent study various soft tissues and water were evaluated. 
In the study by Melhus and Rivard, various tissue-equiv- 
alent materials were compared with water. However, 
it can be concluded that both studies have shown that: 
the phantom material can affect the dose distribution in 
brachytherapy and it depends on the radial distance from 
the source and source's photon energy. 



The results obtained in the present study are geom- 
etry-specific and are valid only for the simulated source 
models. For example the photon spectrum of the model 
6711 125 I source (due to silver (Ag) characteristic X rays) is 
somewhat different from the model 6702 125 I source and 
this phenomenon will affect the dose distribution results 
in various soft tissues. Although some of brachytherapy 
sources such as 192 Ir and 60 Co can be used interchangeably 
in high dose rate (HDR) brachytherapy practice [38,39], 
dose distribution in various soft tissues may differ for 
other lower energy brachytherapy source models, and in 
clinical practice it cannot be sufficient to simply apply the 
results of this study to sources with different energies. 

With interpretation of the relative dose data in Table 6, 
and ignoring a few exceptions, it is evident that the rela- 
tive dose increases with distance from the source model 
for the various soft tissues and water. Thus, it can be con- 
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125 I source 




0.2 0.3 0.4 
Photon energy (MeV) 

169 Yb source 



D 



0.01 0.02 0.03 0.04 

Photon energy (MeV) 

192 Ir source 



0.05 




0.10 0.15 0.20 0.25 
Photon energy (MeV) 



0.4 0.6 
Photon energy (MeV) 



0.2 cm 



1.0 cm 



5.0 cm 



10 cm 



15 cm 



Fig. 3. Photon spectrum at 0.2, 1, 5, 10, and 15 cm distances from the source in 9-component soft tissue. A: 103 Pd source; B: 125 I 
source, C: 169 Yb source, and D: 192 Ir source 



eluded that the effect of tissue composition on dose, or in 
other words the dose difference of soft tissues from the 
dose in 9-component soft tissue is more dominant at dis- 
tal distances from the source. However, when comparing 
the absolute dose values at proximal and distal distances 
from the source models presented in Table 7, one sees the 
absolute dose difference is small, despite the high rela- 
tive difference. Thus, the difference in the compositions 
of various soft tissues on dose has its impact mainly in 
the points close to the source. The effect of variation of the 
relative dose with distance from the source model can be 
related to the energy of photons and photoelectric inter- 
action. This effect is greater for 103 Pd (low energy) source 
model and is less in the case of 192 Ir (high energy) source 
model. Furthermore, when the photons are traversing 
through the material (Table 8), various source models 
have different trends in their spectral changes with dis- 
tance. While for 103 Pd source model the average photon 



energy is increased, for 125 I it is unchanged, for the 169 Yb 
and 192 Ir source models it decreases. The reason for this 
effect may be absorption of a number of low energy pho- 
tons in phantom near the source, for the 103 Pd source. For 
125 I source model, there is a balance between the absorp- 
tion of low energy photons and decreased energy of high- 
er energy photons. For 169 Yb and 192 Ir source models, in- 
teraction of higher energy photons crossing the phantom 
produces a number of low energy photons, which lower 
the average photon energy in the spectrum. 

Figure 3 shows that for all source models there are 
some peaks around more probable energies emitting 
from each radionuclide. It is observed (Figure 3) that in 
the case of 103 Pd, most of points at far distances on the 
curve are above the closer ones. However, for 169 Yb and 
192 Ir source models an inverse trend is observed and for 
the 125 I source model there is a balance between the trends 
for far and near distances from the source. These trends 
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can justify the change of average energy of sources with 
distance as indicated in Table 8. 

As it was plotted in Figure 3, the photon spectrum 
changes with distance from the source at radial direction. 
It can be observed that the number of photons in differ- 
ent ranges of energies decreases with distance. A similar 
pattern is observed for the four sources. This figure can 
provide useful information on the behavior of photon 
spectrum with traverse of the photons inside the phan- 
tom for various brachytherapy source models. While the 
photon spectrum was only illustrated in 9-component 
soft tissue, the same patterns can be expected for other 
phantom materials. The interpretation of this figure can 
also be useful for justification of the dose enhancement 
effect of nanoparticles present in tumor, for tumors at 
various distances from the source [40]. 

In the present study, water as a tissue-equivalent me- 
dia and a number of soft tissues were evaluated based 
on their difference in dose with respect to 9-component 
soft tissue. These media were evaluated as samples to 
illustrate their effect on dose distribution. However, 
there are other soft tissues with different compositions 
and radiation sensitivity. Performing similar evaluation 
on these organs as well as other tissue-equivalent mate- 
rials is a subject for future studies. Furthermore, herein 
a homogeneous and relatively large size of phantom was 
filled with a soft tissue for dose calculations. In a real sit- 
uation, a cross section of the body contains various soft 
tissues and also a bone, while in TPSs the same compo- 
sition is taken for all these tissues. Based on the results 
of the present study, it is predicted that the variation in 
the composition of these tissues will change their doses 
relative to a situation for which the treatment planning 
calculates the dose distribution (a homogeneous case of 
soft tissue or water). A study on a real inhomogeneous 
contour of body, including various soft tissues by Mon- 
te Carlo or 3D dose measurement method (such as gel 
dosimetry) and evaluation of the corresponding error 
compared to a homogeneous soft tissue and the error of 
TPS calculations would be interesting as further works 
in this field. In the present work, the effect of soft tissue 
composition was only studied only for photon energies 
emitted by brachytherapy sources. This effect is import- 
ant in low energy photon ranges, because of dominance 
of photoelectric absorption in low energy range and its 
high dependence on atomic number of the media. 

Conclusions 

In this study, the effect of tissue composition on dose 
distribution for adipose tissue, breast tissue, 4-compo- 
nent soft tissue, brain (grey/ white matter), muscle (skel- 
etal), lung tissue, blood (whole), 9-component soft tissue, 
and water was evaluated for 103 Pd, 125 I, 169 Yb, and 192 Ir 
sources. TPSs, which consider only a homogeneous treat- 
ment volume create errors in dose calculations. Based on 
the results of this study, the error depends on the type of 
soft tissue, brachytherapy source, as well as the distance 
from the source. Furthermore, water which is used wide- 
ly as a tissue equivalent material in radiotherapy dosim- 
etry protocols and TPSs has a dose distribution different 



from other soft tissues when used with photon emitting 
brachytherapy sources. The differences should be taken 
into account in dosimetry protocols and TPSs. 
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